Spinodal decomposition of binary mixtures in uniform shear flow 
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Results are presented for the phase separation process of 
a binary mixture subject to an uniform shear flow quenched 
from a disordered to a homogeneous ordered phase. The ki- 
netics of the process is described in the context of the time- 
dependent Ginzburg-Landau equation with an external veloc- 
ity term. The one-loop approximation is used to study the 
evolution of the model. We show that the structure factor 
obeys a generalized dynamical scaling. The domains grow 
with different typical lengthscales and Ry respectively in 
the flow and in the shear directions. In the scaling regime 
Ry ~ and R^ ~ t"^ , with =5/4 and ay = 1/4. The 
excess viscosity At; after reaching a maximum relaxes to zero 
as 7~^t~^^'^, 7 being the shear rate. Arj and other observables 
exhibit log-time periodic oscillations which can be interpreted 
as due to a growth mechanism where stretching and break-up 
of domains cyclically occur. 
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The kinetics of the growth of ordered phases as a disor- 
dered system is quenched into a multiphase coexistence 
region has been extensively studied in the last years [ 
The main features of the process of phase separa- 
tion are well understood. Typically, domains of the or- 
dered phases grow with the law R{t) ^ where R{t) 
is a measure of the average size of domains. The pair 
correlation function C (r, t) verifies asymptotically a dy- 
namical scaling law according to which it can be written 
as C{r,t) ~ f{r/R), where f{x) is a scaling function. 
In particular, in binary liquids, the existence of various 
regimes characterized by different growth exponents a is 
well established [ |^ . In this letter we study the process 
of phase separation in a binary mixture subject to an 
uniform shear flow. When a shear flow is applied to the 
system, the growing domains are affected by the flow and 
the time evolution is substantially different from that of 
ordinary spinodal decomposition [ ^. The scaling be- 
havior of such a system is not clear. Here we show the 
existence of a scaling theory with different growth ex- 
ponents for the flow and the other directions. For long 
times, in the scaling regime, the observables are modu- 
lated by log-time periodic oscillations which can be re- 
lated to a mechanism of storing and dissipation of elastic 
energy. The behavior of the excess viscosity and other 
rheological indicators reflects this mechanism and is also 



calculated. 

The problem is addressed in the context of the time- 
dependent Ginzburg-Landau equation for a diffusive field 
coupled with an external velocity field [ ^. The binary 
mixture is described by the equilibrium free-energy 
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where ip is the order parameter describing the concentra- 
tion difference between the two components. The param- 
eters b, K are strictly positive in order to ensure stability; 
a < in the ordered phase. The Langevin equation for 
the evolution of the system is 



dip , 



(2) 



where ry is a gaussian stochastic field representing the 
effects of the temperature in the fluid. The fluctuation- 
dissipation theorem requires that 

< ?7(f, t)rj{f', t') >= -2TTV^5{r~ f)S{t ~ t') (3) 

where F is a mobility coefficient, T is the temperature 
of the fluid, and the symbol < ... > denotes the ensem- 
ble average. We consider the simplest shear flow with 
velocity proflle given by 



(4) 



where 7 is the spatially homogeneous shear rate [ ^ and 
Sx is a unit vector in the flow direction. 

In the process of phase separation the initial configu- 
ration of (/3 is a high temperature disordered state and 
the evolution of the system is studied in model (||) with 
a < 0. It is well-know that in this case, also without the 
velocity term, the model (||) cannot be solved exactly [ 

tin this letter we deal with the non-linear term of eq. 
) in the one-loop approximation [ |^J^ . In this approx- 
imation the term ip^ appearing in the derivative 6J^ /dtp 
is linearized as < > ip. It is also called large-n limit. 
Indeed, in the case of a vector field p with n-components 
the term {<p^)<p reduces to < ip'^ > ip in the n — > 00 limit 
[ 1^ . The validity and the limitations of this approxima- 
tion, due to the acquired vectorial character of the order 
parameter, are discussed in literature [ \q. 
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Before presenting our results it is useful to summarize 
the known behaviour of a phase separating mixture un- 
der shear flow. The shear induces strong deformations 
of the bicontinuous pattern appearing after the quench 
[ mii^ ■ When the shear is strong enough stringlikc do- 
mains have been observed to extend macroscopically in 
the direction of the flow [ In experiments a value 
Aa = ax — oiy in the range 0.8 1 for the difference of 
the growth exponents in the flow and in the shear direc- 
tions is measured [ Two dimensional molecular 
dynamic simulations find a slightly smaller value [ p^ . 
We are not aware of any existing theory for the value of 
ax,ay. The shear also induces a peculiar rheological be- 
havior. The break-up of the stretched domains liberates 
an ener gy which gives rise to an increase A77 of the vis- 
cosity [ n4Ulq|. Experiments and simulations show that 
the excess viscosity A77 reaches a maximum aX t — 
and then relaxes to smaller values. The maximum of the 
excess viscosity is expected to occur at a fixed and 
to scale as Ar]{trn) ^ "f^" [ p|,pT[. Simple scahng argu- 
ments predict = 2/3 [ || , but different values have been 
reported [ |l^] . 

We study the time evolution of the structure factor 

C(fc, t) =< ip{k, t)ip{-k, t) > (5) 

where ip(k,t) is the Fourier transform of the field ip{x,t) 
solution of eq. (^. The excess viscosity is defined in 
terms of C(fc, t) by 

Arj{t) = -^-' I kxkyC{k,t) (6) 

J\k\<q {^■^) = 

where g is a phenomenological cutoff. In the one-loop 
approximation the dynamical equation for C{k,t) is: 

- = ~k'[k' + Sit) - l]C(fc, t) + er 

(7) 

where 

/dk 

The parameters F, a, 6, k have been eliminated by a re- 
definition of the time, space and field scales. We solve 
Eq.(^ numerically in two dimensions. A first-order Eu- 
ler scheme is implemented with an adaptive mesh, due to 
the peaked character of the solution. The initial condi- 
tion chosen for the function C(A:,0) is a constant value, 
which corresponds to the disordered state with T = 00. 
The typical evolution of C{k, t) is shown in Fig.l for the 
particular case T — Q and 7 — 0.001. At the beginning 
the function C{k,t) evolves forming a circular volcano 
structure, as usually in the case without shear. This 
is the early-time regime when well-defined domains are 



forming. Then shear-induced anisotropy effects become 
evident in the elliptical shape of C{k,t) and in the pro- 
file of the edge of the volcano, as it can be seen in Fig.l 
at = 0.05. Similar elliptical patterns of C{k,t) are 
usually observed in experiments. The dips in the edge 
of the volcano develop with time until C{k,t) results to 
be separated in two distinct foils, as jt ^ 1. This ex- 
plains the disappearing of the peak corresponding to the 
major axis of the ellipse observed in experiments [ 
During this evolution the support of C{k,t) shrinks to- 
wards the origin with different scales for the shear and 
the flow directions. At later times in each foil of C{k,t) 
two peaks can be distinguished and the relative heights 
of these peaks change in time. In Fig.l at 7/: = 6 the 
peak characterized by \ky\ 3> \kx\ dominates, while the 
other peak with \ky\ ~ \kj:\ prevails successively. The os- 
cillations between the two peaks have been observed to 
continue in time and characterize the steady state. 




FIG. 1. The structure factor at consecutive times for 
7 = 0.001. The kx coordinate is on the horizontal axis and 
assumes positive values on the right of the pictures, while the 
ky is positive towards the upper part of the coordinate plane. 
The support of the function C{k, i) shrinks towards the origin. 
For a better view of C(fc, t), in the last two pictures, we have 
enlarged differently the scales on the kx and ky axes. The 
actual angle between the direction of the foils of C{k,t) and 
the ky axes is 6 = 21° and 9 = 13° in the last two pictures. 

A quantitative measure of the size of domains is 
given by Rx(t) = I/a/ < k^ > where < fc^ 
/ dkkxC{k,t) / J dkC{k,t), and the same for the other 
directions. The evolution of R^jRy is plotted in Fig. 2. 
The growth exponents in the shear and in the flow direc- 
tion are ~ 1/4 and ax — 5/4. The value = 1/4 is 
the same as in models with vectorial conserved order pa- 
rameter without shear; this corresponds to the Lifshitz- 
Slyozov exponent a = 1/3 for scalar fields. A growth 
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exponent ay unaffected by the presence of shear is also 
measured in experiments [ |l^]. We see in Fig. 2 that the 
ampHtudes of R^., Ry oscillate periodically in logarithmic 
time. This behavior can be related to the oscillations of 
the peaks of C(fc,t) previously observed and will be dis- 
cussed later in relation with the behavior of the excess 
viscosity. 
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FIG. 2. The average size of domains in the x and y direc- 
tions as a function of the strain 'yt. The two straight lines 
have slope 5/4 and 1/4. 

In order to study analytically the behavior of the model 
for arbitrary space dimensionality d we resort to a scaling 
ansatz [ |l^ . For the structure factor we then assume 



Cik,t) = Y[R^{t)F \x,T{^t) 



(9) 



for long times, where X is a vector of components Xi — 
kiRi{t), F is a scaling function and the subscript i labels 
the space directions with i = 1 along the flow. We also 
allow an explicit time dependence of the structure fac- 
tor through T{jt)\ notice that since C{k,t) scales as the 
domains volume below the critical temperature, r must 
not introduce any further algebraic time dependence in 
C{k,t). We then argue that F is a periodic function of 
r, as suggested by the oscillations observed numerically 
in the physical observables. Inserting this form of C{k, t) 
into eq. we obtain: 



7X1 ^2 = RiR2^ lidF/dT + J2 
[ i=i 



Rr'R,{F + X,F,y 



RfX^[Y^R-^Xl-l + s]F 



\k=l 



(10) 



where Fi — dF/dXi and a dot means a time derivative. 
Since the l.h.s. of Eq.(|ic|) scales as one has the solu- 
tion R,{t) - 7*'^', T(7t) - log7i, S{t) = 1 - t-^, with 
= 1, ^ (i 2,(i), ai = 5/4, = 1/4 {i = 2,d) 
and /3 = 1/2. In this way we recover the growth expo- 
nents previously found. Actually the exponents found 
numerically are slightly smaller then the predicted pow- 
ers due to logarithmic corrections [ |l^ . 

We now turn to the analysis of the rheological be- 
havior of the mixture and in particular of the excess 
viscosity. The previous theoretical arguments can be 
used to establish the scaling properties of Ary. In- 
serting the form into Eq. (0) we obtain Ar]{t) ^ 



j-^Ri{t)-^R2{t)-\j{T) ^ j-H^/^g{T), where ^(r) 

J X1X2F X,T{t) dX is a periodic function of T{jt). 

Therefore, in the scaling regime, for each value of jt, the 
functions Ar] corresponding to different values of 7 col- 
lapse each on the others if rescaled as Arj — > j-^/^Arf. A 
similar analysis can be done for the normal stress which 

is defined as AA^i = J (2^tt [fey — fc^]C(fc, t) and scales as 

The behavior of Arj at all times, calculated by the nu- 
merical expression of C{k,t), is shown in Fig. 3 for the 
case 7 = 0.001. At] reaches a maximum at jt ~ 3.5, 
then it decreases with the power law t^"^^^ modulated by 
a periodic oscillation in logarithmic time. A comparison 
with Fig. 2 shows that the asymptotic scaling regime 
starts when the excess viscosity reaches its maximum at 
t = trm 3.S found also in experiments [ |ll|. The occur- 
rence of the predicted scaling of A77 with 7 is verified 
numerically with great accuracy for long times. How- 
ever, since tm is at the onset of scaling, an effective ex- 
ponent somewhat larger then 1/2 (i/ ~ 0.6) is measured 
for Ar]{tm), due to preasymptotic corrections. 

The periodic oscillations of A77 are due to the competi- 
tion between the different peaks of C(fe, t). A local maxi- 
mum of A77 occurs for a situation similar to that of Fig.l 
at "ft — 6, when the peak with \ky\ ^ \kx\ dominates 
and the difference between the height of the two peaks is 
maximal. The minima of A77 correspond to the opposite 
situation, as in Fig.l at = 10. The oscillations can be 
explained in this way: The elongation of the domains in 
the flow direction produces an increase of Ar]. Stretched 
domains are characterized by Ry <^ Rx and, therefore, 
are represented by the peak of C{k,t) with \ky \ 3> |fcx|, 
which dominates in this time domain. As time passes, 
however, domains are deformed to such an extent that 
they start to burst dissipating the stored energy. As a 
consequence Aij decreases and more isotropic domains 
are formed. These are characterized by similar values of 
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Rx and Ry and correspond to the other peak of C{k, t). 
This peak starts growing faster than the other until it 
prevails. Later on a minimum of A?/ is observed. Then 
elongation occurs again and this mechanism reproduces 
periodically in time with a characteristic frequency. To 
our knowledge the existence of this periodic behavior has 
never been discussed before 
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FIG. 3. The excess viscosity as a function of the strain 'yt. 
The slope of the straight line is -3/2. 

To conclude, we have studied the phase separation of 
a binary mixture in an uniform shear flow. Dynamical 
scaling holds for this system. Domains grow along the 
flow as Rx{t) ^ i^^^ while in the other directions the 
exponent of the diffusive growth is the same as without 
shear. The difference Aa between the growth exponents 
is 1, a result which is consistent with real experiments. 
The excess viscosity after the maximum relaxes to zero 
as 7~^t~^/^. The amplitudes of physical quantities are 
decorated by oscillation periodic in logarithmic time. It 
would be interesting to study these phenomena in direct 
simulation of the Langevin equation and also to see the 
effects of hydrodynamics on this system. 
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